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ABSTRACT 

Aims. We investigate the dependence of the LMXB population in early-type galaxies on the stellar age. 

Methods. We selected 20 massive nearby early-type galaxies from the Chandra archive occupying relatively narrow range of masses 
and spanning broad range of ages, from 1.6 Gyr to more than 10 Gyrs, with the median value of 6 Gyrs. With the ~ 2000 X-ray point 
sources detected in total, we correlated the specific number of LMXBs in each galaxy with its stellar age and globular cluster content. 
Results. We found a correlation between the LMXB population and the stellar age - older galaxies tend to possess about ss 50% 
more LMXBs (per unit stellar mass) than the younger ones. The interpretation of this dependence is complicated by large scatter and 
a rather strong correlation between the stellar age and the globular cluster content of galaxies in our sample. We present evidence 
suggesting that the more important factor is the evolution of the LMXB population with time. Its effect is further amplified by the 
larger globular content of older galaxies and correspondingly, larger numbers of dynamically formed binaries in them. We also found 
clear evolution of the X-ray luminosity function with age, that younger galaxies have more bright sources and fewer faint sources 
per unit stellar mass. The luminosity function of LMXBs in younger galaxies appears to extend significantly beyond 10 39 erg/s. Such 
bright sources seem to be less frequent in older galaxies. We found that 3 out of x 8 (ultra-) luminous sources are located in globular 
clusters. 

Key words. X-rays: binaries - (Galaxy:) globular clusters: general - Galaxy: stellar content 



1. Introduction 

Low-mass X-ray binaries (LMXBs) are accreting systems con- 
sisting of a low-mass star (<1M ) transferring mass onto a neu- 
tron star or black hole through Roche-lobe overflow. Extensive 
studies of nearby galaxies with Chandra confirmed the long 
suspected fact that their contribution to th e X-ray emission 
in early-type galaxies is substantial (e.g. Ilrwin et a"D l2003t 
iKim & Fabbianoll20 04). Their collective luminosity was found 
to closely follow the near-infrared light and the scaling relation 
of the LM XB population with the stellar mass has been obtained 
dGilfanovN200l . However, a moderate scatter exists in these 
relations, suggesting that the specific number (per unit stellar 
mass) of X-ray binaries is not universally constant among galax- 
ies and secondary correlations may play a role. 

Obviously, one of the main candidates for the second or- 
der correlation is with the age of the stellar population. Unlike 
high-mass X-ray binaries, LMXBs are found both in young and 
old galaxies. Given that the characteristic time scale for the 
stellar evolution of the donor star and for the orbital evolution 
of the binary are both in the Gyrs range, it is not surprising 
that younger and old er galaxies differ in thei r LMXB content. 
For example, recently Kim & Fabbiano] (l2010h reported the evi- 
dence that younger galaxies may have a higher fraction of bright 
sources than the older ones. O n the theoretical sid e, the popu- 
lation synthesis calculations by Fragos^etaTJ d2008l) predict that 
the formation rate of LMXBs steadily decreases with time af- 
ter 1 Gyr. This conclusion seems to have been supported by 
observations - based on the analysis of galaxies dete c ted in 
the extended Chandra Deep Field South iLehmer et al.l d2007l) 
found that for optically faint early-type galaxies (where LMXBs 



dominate the X-ray emission), Lx/Lb increases moderately with 
redshift over z ~ 0.0 - 0.5 range. Other observational facts, 
however, appear to challenge this conclusion. In an SO galaxy 
NGC 5102 , whos e stellar population is younger than 1 Gyr, 
Kraft et al.l (12005b found only two sources brighter than 10 37 
erg/s, whic h is three times less than the predicted number of 
6 LMXBs. Bogdan & Gilfanov (2010) reported similar results 
for two young elliptical galaxies, NGC 3377 and NGC 3585. 
Admittedly, both studies suffered from relatively low statistical 
significance of the results and therefore cannot be con s idered 
as a final argument. Similarly, the result of ILehmer et al.l (120071) 
was not based on a direct determination of the age of the stellar 
environment (which was rather inferred from the redshift) and 
could have been contaminated by other effects (e.g. the rate of 
galaxy mergers). 

It is obvious that more observational effort is needed in order 
to clarify this issue. However, progress in this direction is ham- 
pered by the difficulty in reliable determination of the age of 
stellar populations. In addition, significant fraction of LMXBs 
in elliptical galaxies reside in globular clusters are dynamically 
formed in two-body stellar interactions, rather than having pri- 
mordial origin. In order to investigate the age effects on the pri- 
mordial population of LMXBs, globular cluster sources need 
to be identified and excluded completely from the analysis. To 
this end, clean and reliable lists of globular clusters are needed, 
which are not available for large number of galaxies required for 
a statistically meaningful study. 

By now, Chandra has observed large number of galaxies 
with different morphological types and ages. On the other hand, 
significant progress has been achieved in the accuracy of the age 
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Table 1. The galaxy sample 



Galaxy 


Type 


Distance 


N H 


L K 


MJL K 


r c 


D25 


My 


A'gc 


5 n 






(Mpc) 


(10 20 cm" 2 ) 


(10 10 Lk,o) 


(M /L K , Q ) 


(arcmin) 


(2a, 2b, 8) 


(mag) 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 



N720 


E5 


24.1 


1.54 


16.28 


0.86 


1.20 


4.7', 2.4', 140° 


-21.74 


660 ± 


190" 


1.33 ± 0.38 


N821 


E6 


24.1 


6.39 


9.12 


0.82 


1.66 


2.6', 1.6', 25° 


-21.12 


320 ± 45* 


1.14 ± 0.16 


N1052 


E4 


19.4 


3.07 


8.94 


0.80 


1.12 


3.0', 2.1', 120° 


-21.00 


400 i 


:45 c 


1.59 ± 0.18 


N1380 


SAO 


17.6 


1.31 


12.57 


0.81 


1.32 


4.8', 2.3', 7° 


-21.23 


560 + 


:30 rf 


1.81 ± 0.10 


N1404 


El 


21.0 


1.36 


18.73 


0.85 


0.79 


3.3 ,3.0 , 162.5 


-21.58 


725 ± 


145 e 


1.69 + 0.34 


N3115 


SO 


9.7 


4.32 


9.43 


0.83 


1.07 


7.2 ,2.5 ,40 


-21.13 


630 ± 


150' 


2.22 + 0.53 


N3379 


El 


10.6 


2.75 


7.92 


0.83 


1.17 


£ Al A O / f~ ^1 CO 

5.4 ,4.8 ,67.5 


-19.95 


270 ± 69 s 


2.84 + 0.72 


N3585 


E6 


20.0 


5.58 


18.92 


0.77 


1.20 


4.7', 2.6', 107° 


-21.76 






0.50 ± 0.15* 


N3923 


E4-5 


22.9 


6.21 


29.90 


0.82 


1.66 


5.9', 3.9', 50° 


-22.11 


2494 + 


: 286' 


3.57 ± 0.41 


N4125 


E6 


23.9 


1.84 


23.49 


0.80 


1.95 


5.8', 3.2', 82.5° 


-22.13 






1.30 ±0.50* 


N4278 


El-2 


16.1 


1.77 


7.87 


0.78 


1.15 


4.1', 3.8', 27.5° 


-20.96 


1300 ± 


30O 7 


5.35 ± 1.23 


N4365 


E3 


20.4 


1.62 


20.86 


0.85 


1.66 


6.9', 5.0', 40° 


-22.01 


2511 ± 


1000 J 


3.95 ± 1.57 


N4374 


El 


18.4 


2.60 


24.94 


0.83 


1.70 


6.5', 5.6', 135° 


-22.25 


4301 ± 


1201* 


5.39 + 1.50 


N4382 


SAO 


18.5 


2.52 


27.06 


0.76 


1.82 


7.1', 5.5', 12.5° 


-22.23 


1 1 10 =t 


: 181* 


1.43 ±0.23 


N4472 


E2 


16.3 


1.66 


41.88 


0.85 


3.47 


10.2', 8.3', 155° 


-22.68 


7813 + 


:830* 


6.61 ± 0.70 


N4552 


E0-1 


15.3 


2.57 


10.82 


0.83 


0.98 


5.1', 4.7', 150° 


-21.29 


984 ± 


198* 


2.99 ± 0.60 


N4636 


E0-1 


14.7 


1.81 


13.24 


0.81 


2.95 


6.0', 4.7', 150° 


-21.33 


4200 + 


: 120' 


12.38 ±0.35 


N4649 


E2 


16.8 


2.20 


32.44 


0.85 


2.29 


7.4', 6.0', 105° 


-22.38 


4745 ± 


1099* 


5.32 ± 1.23 


N4697 


E6 


11.7 


2.12 


8.82 


0.77 


2.40 


7.2', 4.7', 70° 


-21.16 


1100± 


400'" 


3.78 ± 1.37 


N5866 


SAO 


15.3 


1.46 


9.47 


0.72 


1.35 


4.7', 1.9', 128° 


-20.93 


400 ± 


100" 


1.69 ±0.42 



(1) - galaxy name. (2) - morphological type. (3) - D istance derived from the surface brightness fluctuation method dTonrv et"aTjl2001l) . (4) 
- Galactic column density (Dickev & Lockman 1990). (5) - Total A^-band luminosity calculated from the total apparent A"s-band magnitude 
from 2MASS Large G alaxy Atlas |jaTOtt*eTaUl2003l) and t he distance adopted in this paper. (6) - A" s -band mass-to-light ratios derived from 
iBell & de Jons (2001), with B - V colors from RC3 catalog dde Vaucouleurs et al.ll 199 lh . (7) - Effective radius from S-band photometry in RC3 
catalog. (8) - D25 region of major diameter (2a), minor diameter (2b), and position angle (9) from RC3 catalog. For NGC 1404, NGC 3379, NGC 
4125, NGC 4278, NGC 4382, and NGC 4552 the position angle is taken from ff s -band image. (9) - Absolute V-band ma gn itude calculated from 
m v from RC3. (10) - Tota l number of globular cl usters. Referen c es - '{K issler-Pati g et al] d!99 6); ' Spit ler et alj d2008T); iForbes et al] (12 001); 
iKissler-Patig et alj dl997l); iForbes et al] dl998h; ^Harris! dl99lh: 4Rhode & Zepfl d2004h ; ISikkema et alj d2006l) ; TForbesl d!996h ; iPeng et al] 
(2008); iDirsch et alj ( 120051) : " iDirschl d!996h ; ICantielloe t al. (2007). (11) - The globular cluster specific frequency calculated from their to- 
tal number and to tal absolute V-band magnitude of the host galaxy, except for the two galaxies for which local values of Sn are used, from 
lHumphrevl ( l2009l) . 



determination techniques and advanced spectroscopical methods 
have been applied to a large number of galaxies. This motivated 
us to undertake a systematic study of the dependence of proper- 
ties of LMXB populations on the stellar age. Among such prop- 
erties we consider the specific (per unit stellar mass) number and 
X-ray luminosity of LMXBs and their luminosity distributions. 
In our analysis we will take into account the possible contami- 
nation by the globular cluster sources to the degree allowed by 
the available globular cluster data. 

The paper is structured as follows: In Sect. [2] we describe 
our selection criteria and the resulting sample. In Sect. [3] we de- 
scribe the X-ray and near-infrared data preparation and analysis. 
In Sect. |4]we discuss average scaling relations for LMXBs and 
their average luminosity functions. Dependence of the LMXB 
numbers and luminosity distribution on the stellar age is dis- 
cussed in Sect. [5] In Sect. [6] we consider the origin of the lumi- 
nous X-ray sources in early-type galaxies, and their dependence 
on the age. Our results are discussed and summarized in Sect. [7] 
and |8] 

2. The sample 

Our goal is to build the largest possible sample covering the 
widest possible range of stellar ages. The size of the sample, 
however, is limited by the content of the Chandra archive and 
by the published age determinations. Our selection criteria are 



the following. Firstly we selected all early-type (E/S0) galaxies 
available in the Chandra archive. We cross-correlated this list 
with publications on the stellar age determinations, leaving only 
galaxies for which reliable age determinations are available (see 
below). From the remaining galaxies we selected only the ones 
located within the distance of ~ 25 Mpc - this ensures a source 
detection sensitivity of better than 5 ■ 10 37 erg/s in less than 150 
ksec of Chandra observation. Then we chose massive systems 
with Lie > 5 ■ 10 10 Lk,o to guarantee the presence of a statistically 
meaningful number of LMXBs (> 20) above the Chandra sensi- 
tivity limit. Finally, we excluded galaxies with ongoing or very 
recent star formation since the stellar content in such galaxies is 
likely to be inhomogeneous. 

In total we selected 20 galaxies with the integrated Zfs-band 
luminosity in the relatively narrow range from ~ 8 • 10 10 to 
4 • 1O h Lk,0. The main properties of these galaxies are listed in 
Table [TJ The Chandra detection sensitivity (Li; m ), which is de- 
fined as the 60% completeness level (Sect. 13.31 ) of LMXBs in 
the study field (Sect. \Jj}, ranges from ~ 4 ■ 10 36 to 10 38 erg/s 
(Table[3]). This ensures that there is statistically meaningful num- 
ber of compact sources in each galaxy. 

2. 1 . The stellar age 

The most accurate and widely used method of the age de- 
termination of elliptical galaxies is the spectroscopic es- 
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timator, which compares observed strength of the ab- 
sorption lines of age-sensitive elements with predictions 
from the simple stellar population (SSP) synthesis mod- 
els. A number of such measurements for diff erent galaxy 
sample s are published i n the l iterature, e.g.. | Trager et all 
(200C j); iKuntschner et ail (120011); ITerlevich & Forbes! (12002k 
Cald well et all (120031): iThomas et al.1 (120051): IDenicolo et al] 
d2005l): ISanchez-Blazauez et al.1 (12006b : lAnnibali et alj (12007b : 
iGallaeher et aTH 2008 ). Several caveats regarding age determi- 
nations are in order. 

What all these measurements have in common is that they 
analyzed the spectrum of the very small central region of the 
galaxy, usually corresponding to r e /8, where r e is the effective 
radius of the galaxy. As LMXBs are rare objects, their num- 
bers inside r e /8 detected in a typical Chandra observation of 
a typical galaxy are by far insufficient for any statistically mean- 
ingful analysis. Quite to the contrary, in order to increase their 
numbers, LMXBs are collected from the region which size is 
comparable with the D25 diameter. Moreover, to avoid compli- 
cations due to the central (weak) AGN, centrally peaked diffuse 
emission and source confusion, the central region of the size of 
~ few arcsec is usually excluded from the X-ray point source 
analysis. Thus, the LMXBs and age measurements are inevitably 
performed in geometrically different, sometimes barely overlap- 
ping regions of the galaxy. Obviously, such analysis requires an 
assumption of homogeneity of the stellar population, which may 
not be fulfilled in all the cases. 

Secondly, there is a (sometimes large) discrepancy between 
different measurements. This discrepancy can have different ori- 
gins. It may be the consequence of the complexity of the task, 
and result from the differences in particular methods used by 
different authors, for example, differences in the procedure of 
correcting for the ionized gas emission, in the choice of absorp- 
tion lines used for fitting or in the libraries of stellar population 
synthesis models. In some cases, these discrepancies may also 
be a consequence of the complexity of the stellar populations in 
the central regions of some galaxies. 

It is known that among different factors, the contamination 
by the gas emission is one of the most important. Therefore for 
galaxies with more than one measurement, we chose the ones 
which treated this issue in a more rigorous way. With this in 
mind, we prior itized the age determ inat ion studies in the follow- 
ing or der: 11 lAnnibali et al.l (j2007l) . 21 ISanchez-Blazauez et alJ 
(20061), 3) IDenicolo et aHfeOOSl) and 4) ITerlevich & Forbes! 
(2002). These 4 papers contain ages for all 20 galaxies in our 
sample. They are summarized in Table [2] where for each galaxy 
we also list the adopted age. 



2.2. The globular cluster content 

To characterize the globular cluster content of a galaxy we use 
the globular cluster specific frequency (Sn), which is c onven- 
tional ly defined by the relation: S N = N GC l0 OA{Mv+15) (iHarrisi 
1199 lb . As this parameter is sensitive to the assumed distance to 
the galaxy, the sensitivity limit and the completeness of the op- 
tical data, we collected the most accurate measurements of the 
total number of GCs {Ngc) in our galaxies, and then computed 
Su with the distances used in this paper. For two galaxies, NGC 
3585 and N GC 4125, Nnr wer e not available and we used the lo- 
cal Sn from Humphrey (2009) as an approximation to its global 
value. All values of Sn are listed in TableQ] 



Table 2. The stellar age measurements. 



Galaxy 


Age 


1 


Age2 


Age3 


Age4 


Adopted age 




(Gyr) 


(Gyr) 


(Gyr) 


(Gyr) 


(Gyr) 


Y N720 








-4-3 


3.4 


->• -4.3 


Y N821 






C O i 1 c 

5.2 ±1.5 


4.0";' 


7.2 


5.2 ±1.5 


N1052 


14.5 ± 


4.2 




9 Q+0-4 




14.5 ± 4.2 


Y N1380 


4.4 ± 


0.7 








4.4 ± 0.7 


Y N1404 










5.9 


5.9 


N3115 






8.4 ± 1.1 


2 6 +0 -° 

-0,6 




8.4 ± 1.1 


N3379 






8.2 ± 1.1 


10.9^ 


9.3 


8.2 ± 1.1 


Y N3585 










3.1 


3.1 


Y N3923 








2 ' 6+ 06 






Y MZL1 9^ 








-'• y -3.0 




• ) - y -3.0 


N4278 






12.5 ± 1.2 




10.7 


12.5 ± 1.2 


N4365 






7.9+1.2 


3.6 + " 




7.9 ± 1.2 


N4374 


9.8 ± 


3.4 


11.3 ± 1.3 


3.8!« 


11.8 


9.8 ±3.4 


Y N4382 








1.6 


1.6 


N4472 






9.6+1.2 




8.5 


9.6 ± 1.2 


Y N4552 


6.0 + 


1.4 


12.4+1.2 




9.6 


6.0 ± 1.4 


N4636 


13.5 + 


3.6 


10.3 ± 1.3 






13.5 ±3.6 


N4649 






16.9 + 2.3 




11.0 


16.9 ±2.3 


N4697 


10.0 + 


1.4 


5.9+ 1.2 




8.2 


10.0 ± 1.4 


Y N5866 










1.8 


1.8 



The st ellar ages are listed according to the prio rity: 1 -lAnnibali et al.1 
j2007l). 2 -ISanchez-Blazauez et al. (2006), 3 - IDenicolo et alJ i2005h . 
4 - ITerlevich & Forbe s ( 2002). All the measurements refer to the central 
r e /8 region. Galaxies marked by "Y" are classified as young galaxies in 
the Sect. 15.21 the remaining galaxies are classified as old. 



3. Data analysis 

3.1. Data preparation and source detection 

Chandra observations of our sample galaxies are listed in 
Table [3] We reduced the data following the standard CIAO 
threads (CIAO version 4.2; CALDB version 4.2.1). We did not 
exclude time intervals for background flares since the benefit 
of the increased exposure time outweighs the increased back- 
ground. The energy range was limited to 0.5-8.0 keV. We made 
exposure maps in this energy range, assuming the single power- 
law model with F = 1 .7 under the galactic absorption for each 
galaxy. To detect point sources we used CIAO task wavde- 
tect w ith the parameters adopted from I Voss & Gilfanovl (2006, 
2007). Thresholds were set to 10~ 6 , yielding averagely 1 false 
detection per 8.2' x 8.2' area (10 6 ACIS pixels). 

For galaxies w ith multiple obser v ations , we corrected their 
offsets following IVoss & Gilfanovl d2007h . using the point 
sources detected in each observation within the 4' radius of the 
telescope axis. The observations were then shifted using CIAO 
task reproject_events, to match the coordinate system of the ref- 
erence (marked with an asterisk in Table 0. The images were 
then combined together and re-analyzed. We performed wavde- 
tect again on the combined images to finalize the point sources 
in each galaxy. To avoid the source crowding problem and the 
bias of the incompleteness of LMXBs in the galaxy center we 
excluded the central a = 5" ellipse region (with the eccentricity 
and position angle following D25). We define the region outside 
the central 5", inside D25 as the study field throughout this pa- 
per. The total number of point sources detected in the study field 
is listed in Table|4] 

To estimate the source counts we applied circular aperture 
centered on the central coordinates (output of wavdetect) of each 
source. We defined the source region as including 85% of the 
local point spread function (PSF) value. The PSF file was ex- 
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Table 3. Chandra observations. 



Galaxy 


Observation ID 


Exposure 


-^-rnin 


Mim 


(1) 


(2) 


(3) 


(4) 


(5) 


N720 


492,7062,7372* ,8448,8449 


138.8 


2.7 


7.2 


N821 


4006,4408,5691,5692,6310, 


212.9 


1.3 


2.8 




6313* ,6314 








N1052 


5910 


59.2 


3.1 


6.3 


N1380 


9526 


41.6 


3.9 


6.1 


N1404 


2942,4174* ,9798,9799 


114.5 


2.4 


11.7 


N3115 


2040,11268,12095* 


153.2 


0.34 


0.70 


N3379 


1587,7073*-7076 


337.0 


0.06 


0.42 


N3585 


2078,9506* 


94.7 


2.3 


4.1 


N3923 


1563,9507* 


102.1 


2.6 


6.3 


N4125 


2071 


64.2 


3.0 


8.9 


N4278 


4741,7077-7081* 


470.8 


0.32 


0.88 


N4365 


2015* ,5921-5924,7224 


195.8 


1.0 


2.5 


N4374 


803,5908*,6131 


115.5 


0.84 


4.9 


N4382 


2016 


39.7 


3.5 


6.3 


N4472 


321*,322,11274 


89.6 


0.58 


5.6 


N4552 


2072 


54.4 


1.3 


4.5 


N4636 


323,324,3926,4415* 


209.8 


0.13 


3.9 


N4649 


785,8182*,8507 


108.0 


2.1 


6.8 


N4697 


784,4727-4730* 


193.0 


0.41 


0.83 


N5866 


2879 


33.7 


2.1 


4.7 



(1) - Galaxy name. (2) - Chandra observation IDs. (3) - The total 
exposure time of Chandra observations. (4) and (5) - The 0.5-8 
keV luminosity of the faintest source detected, and the luminosity 
corresponding to the 60% completeness in the study field (computed 
assuming that the spatial distribution of LMXBs follows the A"„„ 5 -band 
light). The luminosities are in the units of 10 37 erg/s. 



tracted by CIAO task mkpsf from each image, then combined 
together for multiple observations. The background region was 
defined as three times the radius of the source region. For back- 
ground regions having neighboring sources overlapped, we ex- 
cluded their source regions from the background regions. The 
source net counts (with the majority in the source region and mi- 
nority in the background regio n) and errors were then computed 
by the equations (1) and (2) in lVoss & Gilfanovl(l2007l) . To con- 
vert the absorbed source count rates into unabsorbed luminosi- 
ties in 0.5-8 keV, we assumed a power-law spectrum (T = 1.7) 
with galactic absorption. We listed the faintest source detected 
in each galaxy in Tableland the total number of point sources 
above Ly lm in Table [4] 



3.2. The cosmic X-ray background sources 

We estimated the cosmic X-ray background (CXB) sources, us- 
ing the full b and (0.5- 1 ke V) log(AQ-lo g(S ) distribution of CXB 
sources from Georgakakis et al. (2008) and converted the flux to 
the 0.5-8 keV band, assuming a power-law spectrum with a pho- 
ton index of 1 .4. The total number of CXB sources among all 
detected point sources, and point sources above Ln m in the study 
field are listed in Table [4] (The model was corrected by the in- 
completeness function of CXB sources derived in Sect. I3.31 >. In 
most galaxies CXB sources contribute less than 15% of the total 
X-ray population, except NGC 3379 and NGC 4382 where the 
contribution is somewhat higher (25-30%), however the essen- 
tial statistics sustains. 

As well known, the CXB source density is subject to field- 
to-field variations due to the cosmic variance. These variations 
limit the accuracy of the CXB level predictions based on the 
source counts in selected extragalactic fields to ~ 10-30% of the 



predicted CXB value, depending on the solid angle. As the CXB 
contribution to the total number of sources is rather small, these 
uncertainties are relatively unimportant in most of the luminosity 
range. The situation changes in the bright end of the luminosity 
function, where the cosmic variance becomes the major limiting 
factor in our analysis. 

3.3. X-ray incompleteness correction 

The detection sensitivity of point sources varies throughout the 
Chandra images. Various reasons include the non-homogeneous 
level of the diffuse X-ray emission in the galaxy, the deterio- 
ration of the PSF at large off-axis angles, and the non-uniform 
exposure of an image in which observations with different point- 
ings are combined. To calculate the point sou rce detection sen- 
sitivity we used the method and the code from lVoss & Gilfanovl 
(2006), in which the detection method was inverted using the 
local PSF, background and exposure. The incompleteness func- 
tion K(L) is computed as the fraction of pixels weighted by the 
assumed spatial distribution of sources, in which the detection 
sensitivity is better (lower value) than the given luminosity. We 
calculated K(L) separately for the CXB sources and LMXBs, 
since the CXB sources have a flat distribution, while the field 
LMXBs are expected to follow the K rm s -ba nd light (for which 
we used the 2MASS Large Galaxy Atlas data dJarrett et alj2 003) 
- see Sect.[3~4b. 

To estimate the incompleteness-corrected number of 
LMXBs in the study field in each galaxy, we did incompleteness 
correction for the number of all resolved point sources, assum- 
ing K{L) for the field LMXBs, then subtracted the correspond- 
ing number of CXB sources. This procedure is described by the 
equation: 



LMXB 



4nD 



2 dNcxa Kcxb(L) 



dL K L 



dL 



(1) 



where 4nD dNcxs/dL equals to dNcxs/dS, which is the 
log(A0 - log(S) distribution of the CXB sources. We listed the 
total number of LMXBs above L\\ m , CXB subtracted and incom- 
pleteness corrected, in Table [4] 



3.4. Near-infrared data analysis 

We calculated the stellar mass in the study field from near- 
infrared data, using the A"s (2.16 nan) ima ges from the 
2MASS Large Galaxy Atlas dJarrett et al] l2003h provided by 
NASA/IPAC Infrared Science archive. Most images are back- 
ground subtracted, except for NGC 821, for which we obtained 
the background from adjacent regions. We also removed the con- 
tamination of bright fore/background point sources from im- 
ages visually. The integrated, point source- and background- 
subtracted count rate (S ) was converted into calibrated magni- 
tudes with m K (mag) = KMAGZP-2.51og(S), where KMAGZP 
is the zero point magnitude for the ^s-band given in the image 
header. Hence the integrated ^s-band luminosity was calculated, 
and then converted to the stellar mass with the A^-band mass-to- 
light ratios derived from lBell & de Jongld2001l) with B—V colors 
(Table H). We listed the results in TableS 
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Table 4. Statistics of point sources in the study field. 



Galaxy 


jwtotal 


\rtot;il 
' V CXB 


ly'-lim 


JY CXB 


ly'-lim 
JV LMXB 


L K 


M, 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


N720 


79 


5.9 


60 


4.8 


60.8 


14.39 


12.37 


N821 


39 


3.3 


38 


3.1 


36.0 


7.02 


5.76 


N1052 


41 


2.5 


35 


2.2 


35.4 


7.42 


5.93 


N1380 


36 


3.9 


28 


3.4 


28.0 


10.99 


8.90 


N1404 


33 


3.6 


18 


2.5 


17.9 


15.14 


12.87 


N3115 


99 


11.1 


89 


10.2 


82.8 


8.51 


7.06 


N3379 


95 


23.4 


87 


22.9 


65.8 


6.89 


5.72 


N3585 


59 


6.1 


56 


5.8 


53.9 


15.73 


12.11 


N3923 


105 


10.8 


83 


9.8 


89.6 


26.36 


21.61 


N4125 


42 


8.2 


27 


6.6 


24.2 


20.98 


16.78 


N4278 


177 


16.3 


160 


14.9 


154.2 


6.71 


5.24 


N4365 


244 


23.9 


213 


22.5 


201.1 


18.87 


16.04 


N4374 


133 


18.1 


88 


13.9 


91.5 


22.33 


18.53 


N4382 


52 


13.5 


44 


12.3 


33.1 


25.14 


19.11 


N4472 


238 


26.3 


171 


24.0 


171.6 


39.71 


33.76 


N4552 


94 


10.4 


68 


7.6 


70.7 


9.07 


7.52 


N4636 


123 


12.2 


82 


9.1 


83.6 


11.39 


9.22 


N4649 


236 


15.0 


149 


10.7 


168.5 


28.97 


24.62 


N4697 


120 


22.4 


107 


21.4 


85.9 


7.96 


6.13 


N5866 


29 


3.1 


23 


2.7 


21.5 


8.45 


6.09 


Total 


2074 


240.0 


1626 


210.4 


1576.1 


312.0 


255.4 



(1) - Galaxy name. (2) and (4) - Number of all resolved X-ray point 
sources and sources brighter than L\- lm . (3) and (5) - Predicted num- 
ber of CXB sources among (2) and (4). (6) - Total number of LMXBs 
above L\[ m after incompleteness correction and background subtraction. 
(7) and (8) - Total A" s -band luminosity and stellar mass (in units of 10 10 
L K ,o and 10 10 M Q ) in the study field. 




100 
Lx(10 36 erg/s) 

Fig. 1. The observed cumulative distribution of all resolved 
point sources in all galaxies. The distribution is not corrected 
for incompleteness or the contribution of CXB sources. The 
shaded area shows lcr Poissonian uncertainty. The thick solid 
and dotted lines show the predicte d distribution of CXB so urces 
based on the log(AQ - log(5 ) from Georgakak is et al. (2008) and 
iMoretti et all d2003) respectively. 



4. Average XLF and scaling relations for LMXBs 

4.1. XLF of compact sources and the CXB contribution 

The Fig. [TJ shows the combined luminosity distribution of all 
X-ray compact sources detected within the study fields of galax- 
ies along with the predicted distribution of the CXB sources. It 
shows that the CXB contribution is relatively unimportant below 
log(Lx) = 39, where the background AGN account for approxi- 
mately 10% of the observed compact sources. 

Due to quick declining of the LMXB XLF in the log(Lx) ~ 

38.5 - 39 range, accurate account for the CXB contribution be- 
comes crucial at log(Lx) S 39. There is an apparent tail of the 
observed source counts in this luminosity range, which slope is 
similar to the slope of the predicted distribution of CXB sources. 
Its normalization, however, is some what higher than p r edicte d 
by the CXB logfjV) - log(S) from IGeorgakakis et ail J2008). 
Quantitatively, we detected 20 sources above 10 39 erg/s, while 

11.6 background AGN in these fields is predicted based on 
IGeorgakakis et al .1 (120081) . The Poissonian distribution predicts 
a moderately low probability of ~ 1 .6 • 10~ 2 for such a deviation 
solely due to random fluctuations. We also check ed the predic- 
tions o f the CXB log(A0 - log(S) determined by IMoretti et al.1 
(2003) and found that it can fully account for the observed bright 
sources, as shown by the dotted line in Fig.Q] In computing this 
prediction, we used the soft band (0.5-2 keV) coun ts and con- 
verted them to the 0.5 - 8 keV band, as described in Zhan g et alj 
(I2011h . 

Comparing the two predictions, the f ollowing remarks 
shoul d be made. The more recent work of IGeorgakakis et aH 
(2008) is based on a larger sample of sources detected in a larger 
number of Chandra surveys, and its results are in good agree- 
ment with another recent study by iKim et al.l (120071). More over, 
converting the soft band flux from IMoretti et al. r d200l into 
0.5-8 keV highly depends on the assumed spectrum of CXB 
sources, for example, changing the photon index from 1.4 to 
1.9 makes ~ 1.5 times of difference. Thirdly, the numbers of 
bright sources detected outside the D25 region of galaxies tend 
to be in agreement with the prediction of the log(A0 - log(5) 
by IGeorgakaki s et al] (120081) (see be l ow). F or these reasons we 
decided to accept IGeorgakakis et al.l (120081) as our default CXB 
model. 



4.2. Average XLF of LMXBs 

The CXB subtracted and incompleteness corrected cumulative 
X-ray luminosity functions (XLF) of compact X-ray sources in 
each galaxy are plotted in Fig. [2] As low-mass X-ray binaries 
are nearly the only type of compact X-ray sources in early-type 
galaxies capable of emitting at the log(Lx) ^ 36 luminosity 
level, the distributions shown in Fig. [2] can be regarded as lu- 
minosity functions of low-mass X-ray binaries in these galaxies. 
The XLFs have been normalized to unit stellar mass in the study 
field. It is clear that all the XLFs have a similar shape, which 
is broadly consistent with the average XLF of low-ma ss X-ray 
binaries in nearby galaxies obtained by iGilfanovl (120041) (plotted 
with the thick line in the figure). On the other hand, a notable 
scatter of more than a factor of two exists for the normalization, 
which is a manifestation of the scatter in the LMXB - stellar 
mass relation, as discussed below. 

To construct combined XLF of all galaxies with differ- 
ent detection sensi tivity, we followed the method described in 
IZhang et aLld201 lb . The cumulative and differential forms of our 
XLF are plotted in Fig. [3] We fitted the combined XLF with the 
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A 




100 

Lx(10 16 erg/s) 

Fig. 2. Cumulative XLFs of low-mass X-ray binaries in galaxies 
of our sample. The luminosity functions are CXB subtracted and 
incompleteness corrected, and normalized to the stellar mass of 
the host galaxy. They are plotted above corresponding Ln m of 
each galaxy. The sol id line is the av erage XLF of LMXBs in 
nearby galaxies from lGilfanovl d2004l) . 

template introduced in Gilfanov (2004): 



dN 
dLx, 



K2 (L36/Lb,2) " 2 

K3 (iW^cut) - " 3 ■ 




L36 < £b,i 

Lb,l < ^36 < ^b,2 
Lb,2 < ^36 < L cut 

L36 > L cul 



(2) 



where L36 = Lx/ 10 36 erg/s and normalizations ^123 are related 
by: 

K2 = K\ (Lb,i/ib,2) a2 , 

The value of the high luminosity cut-off was fixed at L cut = 5 ■ 
10 4 . We performed Maximum-likelihood fitting to the unbinned 
data. Our best-fit parameters with lcr errors are: a\ = 1.00±0.09, 



a 2 



2.08!^, a 3 



3.67^,1*,, 



582 



+ 107 



i=53+ 4 : « and L b ,2 
The normalization is K\ = 1.1 ± 0.3 per 10 11 M Q . 

The combined XLF obtained in this study is broadly consis- 
tent with the average LMXB XLF obtained by iGilfanovl d2004 
(cf. dotted line in Fig. 0). The XLF of the sources in our sample 
appears to be somewhat flatter in the bright end of log(Lx) > 
38.5, having more luminous sources. As it will become evident 
in the section I5T21 this is related to the large fraction of younger 
galaxies in our sample. Additionally, there is a rather peculiar 
tail of luminous sources above log(Lx) ^ 39, which will be dis- 
cussed in Sect. [6] 



4.3. Scaling relations for LMXBs 

By its definition, the present sample is not designed for detailed 
analysis of scaling relations of LMXBs with stellar mass, as se- 
lected galaxies occupy a rather narrow range of masses, less than 
a factor of 5= 7. The XLF analysis presented above suggests that 
our data is generally consistent with the Nx - M, and Lx - M» 
dependence obtained previously. This is further analyzed below. 

The total numbers and collective luminosities of LMXBs 
more luminous than 5 • 10 37 erg/s (after subtracting the CXB 
contribution and correcting for the incompleteness, Eq.(Q}) for 



x 
_i 

A 




100 

Lx(10 36 erg/s) 

Fig. 3. The combined luminosity distribution of LMXBs in our 
sample in the cumulative (upper panel) and differential (lower 
panel) forms. The distributions are CXB subtracted and incom- 
pleteness corrected, as described in the text. The shaded area 
in the upper panel indicates the lcr Poissonian uncertainty. The 
solid lines show the best-fit model with t wo breaks, the d otted 
lines show the average LMXB XLF from lGilfanovl d2004 . See 
Sect.|6]for the discussion and caveats of the high luminosity tail. 
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stellar mass, 10 10 M Q 

Fig. 4. Relation of the total number (upper panel) and luminosity 
(lower panel) of LMXBs in the study field in each galaxy with 
the stellar mass. 
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Fig. 5. Correlation of the GC specific frequency with the stellar 
age. Filled circles mark galaxies with larger /xlf- 

all the galaxies are plotted against stellar mass in Fig. [4] The 
combined XLF of our sample gives the values of Nx/M* = 54.0 
and L x /M, = 8.2 • 10 39 erg/s per 10 11 M G . The fit to the data 
of individual galaxies shown in the scatter plots in the Fig. [4] 
gives N x /M t = 56.5 and L x /M, = 9.0 ■ 10 39 erg/s per 10 11 M 
These values are by ~ 1 .4 times larger than obtained bv lGilfanovl 
(2004]) - 37.8 and 5.8- 10 39 erg/s. The best-fit relations are shown 
with the solid lines in the plot. The rms deviation of the points 
from best-fit relations are 0.19 dex for Nx and 0.16 dex for Lx- 
The observed scatter is likely caused by the large spread in the 
galaxy ages and (related to it) variation in the contribution of bi- 
naries formed dynamically in globular clusters. These effects are 
further discussed in the next section. 



5. Dependence on the stellar age and the globular 
cluster content 

5. ) . Correlation of LMXBs with age and globular cluster 
content of the host galaxy 

To characterize the number of LMXBs per unit stellar mass in 
each galaxy, we used the following quantity: 



XLF 



N X (L > Zjim) - N C Xfi(L > ^lim) 

M, x f F(L)K LM x B (L)dL 



(3) 



where Nx(L > L/, m ) and Ncxb(L > L/, m ) are the numbers of de- 
tected X-ray sources and predicted CXB sources (Table©, F (L) 
is the best-fit differential XLF and ^lmxb C^O is the incomplete- 
ness function of LMXBs in the given galaxy. The so defined 
fxLF is the specific (per unit stellar mass) XLF normalization 
computed from the number of resolved LMXBs above L]j m , as- 
suming the average XLF shape of the resolved sources modified 
by the point source detection incompleteness. For the average 
XLF we used the best-fit model from Sect. [4] /xlf equals to 1 
with a mean number of 5.4 LMXBs with luminosity Lx > 5 • 10 37 
erg/s per 10 10 M Q . 

Obviously, /xlf has the advantage that all detected X-ray 
sources above 0.6 incompleteness level are involved in the cal- 
culation. The disadvantage is that it relies on the assumption that 
the XLF shape does not change from galaxy to galaxy. To inves- 
tigate the importance of the latter, we computed /xlf for several 



galaxies using the average luminosity functions for young and 
old galaxies (see the next section) and found that they differ by 
no more than ~ 30%. This is insignificant given the scatter of 
the points and the amplitude of the correlations found below. 

The study is complicated by existence of a strong correlation 
between the stellar age and the specific frequency of globular 
clusters, as illustrated by the Fig. [5] Although we did not find 
any report about such correlation in the literature, its existence is 
not surprising as more GCs are expected in older galaxies, since 
more mas sive clusters a re formed in larger star bursts at larger 
redshifts (Bastian 2008). A significant fraction of these will sur- 
vive through the following evolution of the galaxy (Fall & Zhang 
2001). Interestingly, we do not see any notable correlation be- 
tween the Sn and st ellar mass o f the ga laxy, for example similar 
to the one found by iPeng et alj (12008). This may be a result of 
large scatter in these correlations. A detailed investigation of the 
behavior of the specific GC frequency and its dependence on the 
environment and formation history of the host galaxy is beyond 
the scope of this paper. However, presence of the SiM-age corre- 
lation (at least in our sample) significantly complicates the anal- 
ysis of the dependence of the LMXB population on the stellar 
age and the globular cluster content of the galaxy. 

We plot /xlf versus the stellar age and the GC specific fre- 
quency in Fig. [6] Both plots show moderate trends that the /xlf 
increases with the stellar age of the host galaxy and its globu- 
lar cluster content. Despite rather large scatter of the points, the 
Spearman rank-order correlation test gives the null hypothesis 
probability of p = 0.009 for the age- and somewhat large value 
of p — 0.017 for the S n -dependence. These numbers indicate a 
moderately significant, w 2.5 - 3<x, correlation. 

The correlation of f xi f with S Nis similar to the one found in 
previ ous studies (e.g. Humphrey & Buote 2008| iBoroson et al] 
2011). The presence of such correlation was interpreted as an 
evidence suggesting that significant part (if not all) of the LMXB 
population, including the field sources, was formed in globular 
clusters and subsequently expelled into the field. The existence 
of the equally strong correlation of the specific LMXB frequency 
with the age of the stellar population suggests that this interpre- 
tation is not unique and the more important correlation may be 
with the stellar age. 

In order to investigate this further, we fit the data with a two- 
parameter linear model in the form /xlf = axt + bxSN + c. 
Using x 1 minimization we found the following values of best-fit 
parameters: a = 0.040 + 0.008, b = 0.049 + 0.014, c = 0.399 + 
0.049, with a very large value of^ 2 = 144.6 for 17 d.o.f. With 
these best-fit values the contribution of the age term is about two 
times larger that the contribution of the term, suggesting that 
the more important parameter is the age, rather than the globular 
cluster content of the host galaxy. However, because of the rather 
large dispersion of the points and correspondingly large value of 
X 1 any firm conclusion is premature. 

The average values of /xlf for galaxies younger and older 
than 6 Gyrs are 0.74 + 0.07 and 1.04 + 0.06 respectively. The 
statistical significance of the difference between these two values 
is w 3.2cr, in agreement with the Spearman test results. The total 
number of LMXBs per unit stellar mass above 5 • 10 37 erg/s in 
the young sample is 4.18 + 0.27, which is -2/3 of that in the old 
(6.27 + 0.26). The prediction from the average XLF is 5.4. 

5.2. XLFs of young and old galaxies 

In order to investigate the age dependence of the LMXB XLFs, 
we divided galaxies into young and old groups using their me- 
dian age (6 Gyrs) as a boundary. Each group contains 10 galaxies 
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Fig. 6. Correlation of the XLF normalization, /xlf (Eq. HJ, with the stellar age (left) and the specific frequency of globular clusters 
(right). Circles mark galaxies with Sn > 2.5 (left panel) and rectangles mark galaxies with the stellar age >6 Gyr (right panel). 



which are marked correspondingly in Table[2] The study regions 
in young and old galaxies cover a total solid angle of 125.9 and 
251.6 arcmin 2 respectively, with a total stellar mass of 1.24 and 
1.32-1O 12 M . 

We constructed combined XLF of each group and have them 
plotted in Fig. [7j In general older galaxies have deeper Chandra 
observations which have reached a sensitivity of ~ 5 • 10 36 erg/s, 
while the young group has a sensitivity of ~ 3 • 10 37 erg/s. The 
overall shape of the XLF for young galaxies is flatter than that 
of the old ones. This be havior is in agreement with findings of 
iKim & Fabbianol(l2010l) . 

Similarly, we use the median value of the Sn distribution 
(Sn = 2.5) to divide galaxies into GC-rich and GC-poor sub- 
groups. The resulting XLFs are shown in Fig. [8] As expected 
(from the age dependence of the XLF and the general correlation 
between age and Sn), the combined XLF of GC-rich galaxies is 
steeper than that of the GC-poor ones. However, unlike XLFs of 
young and old galaxies, they appear to be rather similar in the 
bright end log(L x ) > 38.5 - 39. 



6. The nature of the (ultra-) luminous X-ray sources 

After subtracting the contribution of background sources, a 
rather peculiar tail of luminous sources with \og(Lx) ^ 39 re- 
mains (Fig. [3). Although the CXB contribution is insignificant 
throughout the most of the luminosity range, its correct account 
and subtraction is critical for establishing the nature of bright 
sources. In total, we detected 20 compact sources with luminos- 
ity exceeding 10 39 erg/s. Subtracting f rom this number the 11.6 
background AGN predicted from the iGeorgakakis et al.l ([2008) 
CXB log(A/) - log(S), we obtain that ~ 8 - 9 sources should 
be associated with galaxies from our sample. Similarly, among 
sources with Lx > 2 ■ 10 39 erg/s (7 sources observed, 4.2 back- 
ground AGN predicted) ~ 2 - 3 sources are expected to be as- 
sociated with galaxies. Assuming Poissonian distribution, these 
numbers correspond to moderately low probabilities of being a 
result of pure statistical fluctuations: ~ 1.6- 10~ 2 and 0.13 for the 
two luminosity ranges, thus suggesting that the luminous sources 
may indeed be X-ray binaries. 

On the other hand, it is well known that the CXB source 
counts produce somewhat different results in different sky fields 
due to the cosmic variance. Moreover, the slope of the bright tail 




100 1000 
bc(10 36 erg/s) 

Fig. 7. The X-ray luminosity functions of LMXBs in the young 
and old galaxies in cumulative (upper panel) and differential 
(lower panel) forms. The data for old galaxies (red in the color 
version of the plot) is marked by circles in the lower panel and is 
surrounded by the shaded area showing the 1 <x Poissonian uncer- 
tainty in the upper panel. Statistical uncertainty for young galax- 
ies has comparable amplitude. 



of the luminosity distribution of all compact sources in Fig.Q]is 
similar to the slope of the CXB log(TV) - log(S). It is therefore 
possible that tail of bright sources in Fig.[3]is due to background 
AGN unaccounted for due to cosmic variance. The amplitude of 
CXB source density variations depends primarily on the consid- 
ered angular scales, decreasing as the solid angle increases. As 
we combined data of 20 galaxies distributed over extragalactic 
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Table 5. Sources brighter than 10 39 erg/s. 



Galaxy 


RA(J2000) 


DEC(J2000) 


Luminosity 


D25 


opt. 


Young Galaxies (< 6 Gyrs) 


IN /ZU 


i n i . .n i n 
+U1. j j.kjLAZ 


-13.44. Ly.DJ 


1 K 


+ U.Uo 


U.1U 




IN 1 joV 


+\JJ. j0.2O.j0 


-j4. J0.jo.y0 


i no 


_i_ n no 

± u.uv 


U. / j 








■Jt.j7. 10.U7 




+ 04 


47 




N1404 


+03:38:51.99 


-35:35:59.93 


1.14 


±0.07 


0.20 






+03:38:54.78 


-35:35:00.96 


1.21 


±0.07 


0.57 




N3923 


+ 11:50:58.65 


-28:49:13.16 


1.31 


±0.08 


0.38 






+ 11:51:09.54 


-28:48:00.67 


2.98 


±0.12 


0.68 






+ 11:51:06.22 


-28:46:49.91 


3.50 


±0.13 


0.65 


Q 


N4125 


+ 12:08:07.46 


+65:10:28.61 


7.41 


±0.23 


0.05 




N4382 


+12:25:20.32 


+ 18:13:01.41 


1.12 


±0.09 


0.58 






+12:25:17.17 


+ 18:13:46.52 


3.76 


±0.18 


0.93 


Q 


N4552 


+ 12:35:45.77 


+ 12:33:02.46 


1.14 


±0.06 


0.61 


GC 




+ 12:35:41.22 


+ 12:34:51.43 


1.18 


±0.06 


0.62 


GC 


Old Galaxies (> 6 Gyrs) 


N3379 


+ 10:47:50.01 


+ 12:34:56.77 


2.14 


±0.03 


0.04 




N4365 


+ 12:24:26.36 


+07:16:53.55 


1.53 


±0.05 


0.71 




N4374 


+ 12:25:11.92 


+ 12:51:53.81 


10.38 


±0.16 


0.74 


Q 


N4472 


+12:29:41.00 


+07:57:44.46 


1.96 


±0.08 


0.62 






+ 12:29:34.46 


+07:58:51.63 


1.33 


±0.06 


0.79 


G 




+ 12:29:42.33 


+08:00:07.96 


1.02 


±0.05 


0.25 


GC 


N4649 


+ 12:43:46.90 


+ 11:32:34.19 


1.54 


±0.06 


0.48 





Columns are the host galaxy name, coordinates, luminosity in units of 10 39 erg/s, offset from the center in units of the D25 radius, optical 
counterpart from NED (Q: Quasi-stellar object, GC: globular cluster, G: galaxy). 



sky, we do not expect cosmic variance to be particularly strong 
in the combined LMXB XLFs shown in Figs. [3] and [7] In order 
to investigate this further, we checked the source numbers out- 
side the D25 but within ~ 10' x 10' around the aim point of the 
Chandra observation. In total 38 sources with log(Lx) > 39 were 
detecte d, while the CXB log(JV) - log(S) of iGeorgakakis et all 
(2008) predicts 47.8 background AGN. These two numbers are 
consistent within ~ 1 .4cr, which suggests that the effect of cos- 
mic variance is not very strong. Moreover, the observed local 
CXB source density is ~ 20% lower than predicted, i.e. in com- 
puting the LMXB XLF shown in Fig. [3] the CXB contribution 
may be somewhat over-subtracted. 

Furthermore, we checked the luminous, Lx > 10 39 erg/s, 
sources individually. In particular, we looked up their identifi- 
cations in the NASA Extragalactic Database and checked their 
location relative to the D25 radius of the host galaxy. The re- 
sults are presented in Table[5] In the table, we list young and old 
galaxies separately, as before using the median age of 6 Gyrs 
as the boundary. Among the 20 luminous sources, 13 are de- 
tected in the young galaxies and 7 in the old ones, with the 
CXB source predictions being 5.6 and 6.0 respectively. Taken at 
the face value, these numbers suggest that all luminous sources 
are located in young galaxies. On the other hand, the luminous 
sources are evenly distributed between GC-rich (12 sources de- 
tected, 7.2 background AGN expected) and GC-poor galaxies (8 
sources detected, 4.4 background AGN expected). 

According to the NED, four of the luminous sources are 
confirmed background AGN, including the brightest sources 
with L x ~ 10 40 erg/s located in NGC 4374. As expected, 
this number is significantly smaller than the predicted num- 
ber of ~ 12 AGN, due to incompleteness of the NED. 
Interestingly, three luminous sources are located in globular 
clusters, one of which is in the old sample. All three sources 
have luminosities ~ 10 39 erg/s which makes them plausible 
black hole candidates. The source in NGC 4 472 has been 
previously reported by Maccarone et al ] (120071) . One source 



(+12:35:45.77 +12:33:02.46) in NGC 4552 was reported by 
IXu et al] (12005b . The other one (+12:35:41.22,+12:34:51.43) in 
NGC 4552, to our knowledge, is a new identification. 

The above analysis suggests that luminosity function of com- 
pact sources in young sample extends significantly beyond 10 39 
erg/s, with the specific frequency of luminous sources » 6 + 3 
sources per 10 12 M Q . Although the numbers of detected sources 
and uncertainties of the CXB subtraction do not allow us to 
constrain its shape, it appears to be a smooth extension of the 
LMXB XLF observed at lower luminosities. The luminosities of 
these sources, up to ~ 3 - 7 • 10 39 erg/s, are compatible with 
or exceed only slightly the Eddington limit for a ~ 10 M black 
hole. We therefore believe that most likely these sources are stel- 
lar mass black holes accreting from a low or intermediate mass 
companion. It was noted that in the case of an evolved com- 
panion such systems m ay be transient sources dKing et al]|1997t 
iPiro & Bildsten 2002). In order to check this possibility we in- 
vestigated variability of those luminous sources for which mul- 
tiple observations were available. Although bright sources show 
significant variability, up to a factor of ~ 5, we did not find any 
evidence for transient behavior. 



7. Discussion 

One of the two difficulties of this study is the uncertainties and 
limitations of the stellar age determinations. These uncertain- 
ties are the reason for the discrepancy between the values of the 
stellar age determined by different authors. Moreover, age mea- 
surements used in this paper refer to the central few arcsec of the 
galaxy (r e /8). Although the stellar content in elliptical galaxies 
is expected to be more uniform than in late-type galaxies, inho- 
mogeneities are known to exist and in some cases are known to 
be rather large. In order to increase the number of sources, the 
X-ray measurements are carried out over a much larger region 
of the galaxy than where the age was determined. This may give 
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Fig. 8. The X-ray luminosity functions of LMXBs in galaxies 
with high and low globular cluster content Sn (see text for de- 
tails), in cumulative (upper panel) and differential (lower panel) 
forms. The data for high-SN galaxies (red in the color version 
of the plot) is marked by squares in the lower panel and is 
surrounded by the shaded area showing the Icr Poissonian un- 
certainty in the upper panel. Statistical uncertainty for 1ow-5n 
galaxies has comparable amplitude. 



rise to significant noise in the LMXB - age correlation produced 
in this paper. 

To investigate this further, we collected from the literature 
the results of long-slit spectroscopy of galaxies in our sample. 
We have found that some galaxies may have more complicated 
structure of the stellar population than a single-aged population, 
and that the true stellar ages may be different from the values 
used in this paper. We discuss these results below galaxy by 
galaxy. NGC 720 was found to be formed by a merger of an 
old (13-5 Gyr) small-scale s pheroid and a young er (5-2.5 Gyr) 
large-scale disk component dRembold et al.ll2005l) . Thus the av- 
erage stellar age of NGC 720 should be older than 3 Gyr value 
produced by the simple spectroscopic estimators and used in this 
paper. In NGC 821 the central stellar population was found to be 
young (~4 Gyr), while the ag e increases outwar ds and reaches 
12 Gyr at one effective radius dProctor et al.ll2005t) . Thus the av- 
erage age is likely to be older than 5 Gyr. NGC 4125 was found 
to experience a recent dissipational merger event, th us young 
stella r population quite possibly exists in the galaxy (|Pu et alJ 
120101) . And NGC 4365 is very much likely to be older than 7.9 
Gyr, that the decoupled core and the main body of the galax y 
was found to have the same age of ~ 14 Gyr dDavies et alfeoOlh . 
With these numbers, in the left panel of Fig. [6] the data points 
for NGC 720, NGC 821 and NGC 4365 may shift to the right, 
while NGC 4125 may shift to the left. Other galaxies either do 
not have obvious stellar age gradients, or do not have such de- 
tailed individual studies published. For the statistical point of 
view, a study based on the Sloan Digital Sky Survey data for 



~50000 nearby galaxies found that, no age gradient is present 
for massive early-type galaxies (> 1O 1O M ) with a central age 
older than 6 Gyr, while for the ones with central age younger 
than 6 Gyr age gradients can be as big as V age ~ 0.4 (2.5 times 
difference) (To rtora et al.l20l6T) . This result suggests that ages of 
old galaxies in our sample should be on average sufficiently re- 
liable whereas an age uncertainty of up to 0.4 dex is possible for 
galaxies with the central age younger than 6 Gyrs. 

The second difficulty is the correlation between the stellar 
age and the globular cluster content. We did not address the uni- 
versality of this correlation, but the fact that it is rather strong 
in our sample (FigJSJl complicates the separation of the effects 
of these two factors on the LMXB population (Fig|6]l. Results 
of the linear regression analysis suggest that the more important 
factor is the stellar age, with the role of the specific frequency of 
globular clusters being less significant. However, the large scat- 
ter of the points and correspondingly large values of x 1 do not 
allow us to make a statistically sound conclusion. 

The X-ray binaries dynamically formed in globular clusters 
can reside in their parent clusters or can be expelled into the 
field, either due to kicks and/or tidal interactions or in the course 
of destruction of the parent cluster. The former could in principle 
be identified and removed from the sample, provided that clean 
and complete globular cluster lists are available. This however is 
not possible, due to limitations of the globular cluster data, sig- 
nificant part of which was obtained by the ground based facilities 
and suffers from incompleteness and foreground star contamina- 
tion. Although for a number of galaxies HST data is available, it 
also has limited completeness. Besides, in the luminosity range 
of interest X-ray source in globular cluster do not account for 
more than - 30 - 40% of the LMXB population (Zhan g et all 
which is insufficient to explain the amplitude of the ob- 
served dependence (right panel in Fig[6]i. The contribution of 
binaries, dynamically formed in globular clusters and expelled 
into the field is unclear. Based on the correlations of the specific 
frequency of LMXBs with that of globular clusters it has been ar- 
gued that entire population of L MXBs has been f ormed dynam- 
ically in globular clusters (e.g. IWhiteetal.1 [2002b . Our results 
however suggest that the more important factor is the evolution 
of the LMXB population with time. A part of the correlation of 
the specific frequency of LMXBs with that of globular clusters is 
real, due to the contribution of sources formed dynamically (and 
mostly retained) in globular cluster. Another part is however an 
artifact of the dependence of the globular cluster content on the 
age of the galaxy. 

This scenario is further supported by the comparison of 
the luminos i ty fun ction of young and old galaxies (FigJTJ. 
IZhang et"a"D d201 lb has shown that the XLF of globular clus- 
ter sources in the log(Lx) ~ 37 luminosity range is flatter than 
that of the field sources (e.g. see their Fig. 7). Therefore, as older 
galaxies have larger globular cluster content, i.e. larger fraction 
of dynamically formed binaries, their XLF should be expected to 
be flatter than the XLF of the young sample. This prediction is at 
odds with our observations - the XLF of old galaxies is steeper 
that XLF of the young sample. This suggests that the effects of 
evolution of the LMXB population are stronger. 

8. Summary 

The main goal of this paper is to study dependence of the popu- 
lation of low-mass X-ray binaries on the stellar age. To this end, 
we have collected 20 nearby early-type galaxies which were ob- 
served by Chandra to sufficient depth and had the stellar age 
measured. 
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1. We found that older galaxies tend to host more LMXBs per 
unit stellar mass that the younger ones (FigO. The corre- 
lation has large scatter, with the points occupying the dy- 
namical range of a factor of 4. When averaged over young 
(f < 6 Gyrs) and old (f > 6 Gyrs) sub-samples, the specific 
frequency of LMXBs with Lx > 5 • 10 37 erg/s varies from 
4.18 ± 0.27 to 6.27 + 0.26 per 10 10 M Q . Interpretation of this 
dependence is complicated by the rather strong correlation 
between the globular cluster content of the galaxy and its 
stellar age. We presented evidence suggesting that the more 
important factor is intrinsic evolution of the populations of 
LMXBs with time. Its effect is further enhanced by the larger 
globular cluster content of older galaxies resulting in larger 
numbers of dynamically formed binaries in them. 

2. There is clear evolution of the X-ray luminosity function 
with age - the XLF of older galaxies is steeper in the en- 
tire studied luminosity range, log(Lx) ^ 37.5, than the XLF 
of younger galaxies (FigfTJ- 

3. Young galaxies host a significant population of (ultra-) lu- 
minous X-ray sources with luminosity exceeding 10 39 erg/s. 
We estimate their specific frequency of » 6 ± 3 sources per 
10 12 M in the young sub-sample. Such sources are signifi- 
cantly less frequent in the old sub-sample 1 source against 
« 7), with the 90% upper limit of « 2.9 sources per 10 12 M 

4. As a byproduct of this study we found 3 black hole can- 
didates in globular clusters, of which two was previously 
known and one is identified for the first time (Table|5]). 
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